Persistent currents in a circular array of Bose-Einstein condensates 
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A ring-shaped array of Bose-Einstein condensed atomic gases can display circular currents if 
the relative phase of neighboring condensates becomes locked to certain values. It is shown that, 
irrespective of the mechanism responsible for generating these states, only a restricted set of currents 
are stable, depending on the number of condensates, on the interaction and tunneling energies, 
and on the total number of particles. Different instabilities due to quasiparticle excitations are 
characterized and possible experimental setups for testing the stability prediction are also discussed. 
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I. INTRODUCTION 

The existence of flow without dissipation is a signa- 
ture of superfluidity and superconductivity. This phe- 
nomenon has been studied in various condensed-matter 
systems, most notably in helium and metallic super- 
conductors. The experimental achievement of Bose- 
Einstein (BEC) condensation in alkali-metal atomic gases 
opened the possibility of investigating the superfluid 
properties of dilute BEC gases trapped in various ex- 
ternal potentials generated by magnetic or laser fields. 
The geometry of the trapping external potentials can be 
engineered into various shapes: from almost spherically 
symmetric to highly anisotropic, from double wells to lat- 
tices. 

We study, in the present paper, the stability of cur- 
rents for a system of M condensates in a topologically 
restricted configuration, with tunneling allowed only be- 
tween neighboring condensates. This type of trap can be 
created, for instance, for not too large M's, by using an 
experimental setup similar to that of the MIT interfer- 
ence experiment [2| but with a spatial shape of the blue- 
detuned laser beam, which would ensure the separation 
of the gas into three or more pieces. Another possibil- 
ity involves the use of two-dimensional optical lattices in 
magnetic traps of the same type as those recently em- 
ployed to study the superfluid-Mott insulator transition 
3] . Overlapping a highly repulsive optical potential (for 
example, generated by a blue-detuned cylindrical laser 
beam) in the middle of the magnetic trap results in the 
creation of a Mexican-hat potential for the atoms already 
confined in the lattice. 

A simple application of the concept of spontaneously 
broken gauge symmetry to this system has as a conse- 
quence the appearance of well-defined phase relations be- 
tween consecutive condensates and the flow of a circular 
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current around the ring. Nowadays there is a variety 
of experimental methods for inducing currents into an 
already formed condensate: tilting the lattice in a gravi- 
tational field Q or accelerating it displacing the en- 
veloping external potential created by magnetic trapping 
@, phase imprinting 0, coherent conversion between 
two hyperfine states (one serving as the pinning poten- 
tial) y| , or stirring the condensate with a laser beam [jj . 

In this paper we show that, no matter what mechanism 
is chosen to attempt to generate these circular flow states, 
some of them are in fact unstable, either dynamically or 
thermodynamically. The states that are stable and pro- 
duce persistent currents are identified as those with cir- 
culation below a certain critical value, which in general 
depends on the parameters related to on-site interaction 
such as density and scattering length, on the tunneling 
rate, and on the number M of condensates. The maxi- 
mum critical circulation (vorticity) of these states is M/4. 

In the following we will briefly discuss the relevance 
of our study for the two processes that involve at dif- 
ferent points the concept of broken phase symmetry: the 
Kibble-Zurek mechanism and the non-adiabatic quantum 
phase transitions. 

The Kibble-Zurek mechanism ^(J predicts the appear- 
ance of topological defects in systems undergoing a rapid 
quench through a continuous phase transition. For the 
case of dilute alkali atomic gases it has been already ar- 
gued that during a fast quench vortex lines will be formed 
11]. It is likely however that a Kibble-Zurek experiment 
with atomic gases will not be realized in the way origi- 
nally imagined, but rather by crafting traps simulating 
a closed array of condensates which can tunnel from one 
site to another, in a way similar to the experiments done 
using Josephson junctions 01- This is precisely the sys- 
tem analyzed in this paper. A fast evaporative cooling 
of a gas of bosons in this type of trap would result in the 
formation of domains with broken phase symmetry in the 
wells of the potential (due to higher particle densities). 
The overlap of the order parameters leads to transfer of 
atoms between neighboring wells. But in the end, the 
probabilities assigned to different outcome states by any 
microscopic analysis (time-dependent Ginzburg-Landau 
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GL, kinetic theory, etc.) of the non-equilibrium problem 
have to be supplemented with the constrain that only 
the stable currents will survive - indeed, as shown in this 
paper, modes that otherwise would pass unsuppressed 
into the order parameter can decay by the emission of 
quasiparticles. 

A similar conclusion can be drawn for the second pro- 
cess mentioned above, namely quantum phase transi- 
tions. The superfluid - Mott insulator transition has been 
achieved and studied experimentally intensively in recent 
times in two-dimensional optical lattices Q; the same 
ideas can be applied for an array of the type described 
in this paper. In this case crossing the quantum critical 
point results in a coherent superposition of states with 
different circulations 0|. The mechanism that breaks 
the phase symmetry is then any decoherence process: a 
measurement, particle losses, etc. But, roughly, as shown 
below, only half of the macroscopically occupied states 
that result in this way will be persistent currents. 



II. A SIMPLE MODEL 

Let us consider a simple model that captures the essen- 
tial stability features of a macroscopically occupied mode 
in the topologically constrained configuration described 
above. 

We consider M identical small condensates in contact 
with each other and a cylindrical system of coordinates 
(r,9,z). The total number of atoms is N. The centers 
of the condensates are positioned at R\ = (R, 9\,z = 0), 
where 9\ = 2-kX/M and A runs from to M — 1. 

The Hamiltonian of the system is 



H = J drft(r) 



—A 
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where g — Air h 2 /m and a is the scattering length. 

The external potential V(r) is crafted to be high 
enough around the origin, so that the atoms cannot pene- 
trate there; it also has a number of M minima at {R, 9\), 
around which the condensed atoms tend to localize. The 
derealization effect comes from the possibility of tun- 
neling between nearby wells. We also assume that each 
condensate has a small enough number of particles so 
that they are not in the Thomas-Fermi regime (for wells 
of dimension rs 1/im and the scattering length of Na or 
Rb, the number of particles in each condensate can be at 
most in the hundreds). In this situation, the wavefunc- 
tions of each condensate depend weakly on the number 
of atoms in the well 0], and one can apply a M-mode 
approximation for the field operator 



M-l 



W)= ^-R x )a{\), 



(2) 



A=0 



where <p is a solution of the Schrodinger equation for each 



well; the Hamiltonian takes the Bose-Hubbard form 



M-l 



H = -t [a t (A)a(A + l) + a t (A + l)a(A)] 



w 



A=0 
M-l 



« t (A)a t (A)a(A)a(A). 

A=0 



(3) 



Here the constant terms are omitted; the tunneling ma- 
trix element is given by 



t 



V{r) 



dr(j)*{r-R x ) 
and the on-site energy is 

w = g f dr\(j)(r)\ 4 



-Rah 



(4) 



(5) 



The boundary conditions are periodic by construction, 
d(A) = d(A + M). The M-mode approximation and the 
relation can be seen also as an expansion of the field 
operator in the Wannier basis or a tight-binding approx- 
imation. 

The Hamiltonian can be diagonalized in the Bogoli- 
ubov formalism. We start first with identifying the solu- 
tions of the Gross-Pitaevskii equation. This can be done 
by seeing the a's as field operators and A as a discrete 
coordinate. Then the first term of the Hamiltonian © is 
the equivalent of the kinetic energy and the second term 
is the interaction energy. To find the Gross-Pitaevskii 
state we expand the field operator into M modes 



a ( A ) = Y Xk ^ bk ' 

k 



(6) 



where the modes Xfe(A) are periodic Xfc(A) = Xfc(A + M), 
and k takes M positive and negative integer values in 
the interval (— M/2, M/2]. We minimize the mean-field 
energy on a state (Nl) -1 / 2 ^^ \ 0) with the restriction that 
the total number of particles is constant 
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M 



(lf)-M^(at(A)fl(A)) 



A=0 



0. 
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This gives the one-dimensional lattice Gross-Pitaevskii 
equation 

-t [ Xq (X + 1) + x 9 (A - l)] + Nw\ Xq (\)\ 2 x q W = 

(8) 

The analogy with the usual continuous Gross-Pitaevskii 
equation is even more transparent if one notices that if 
x(A) does not vary much from site to site, one can define 



^ X (A) ~ X (A + 1) + *(A — I) — 2 X (A). 



(9) 



We can also write the Hamiltonian (|3J in the momentum 
basis, defined by 



<M 



M-l 

£< 

A=0 



-i(2n/M)k\ 



a(A), 



(10) 
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bl = 



/M)feA-f 



«t(A). 



(11) 



A=0 



The result is, with the indices fc, k', and ^ taking all 
integer values in the interval (— M/2, M/2], 

H = -2i]Tcos f^k) blh + E " b l + i" b lJk>b k , 



k,k'l 



which is formally the Hamiltonian for a uniform system of 
free bosons with kinetic energy — 2t cos (2nk/M) instead 
of the usual h 2 k 2 /2m. 

Equation (JSJ has solutions 



X,(A) = 



*M 



2n N 

u„ = — 2tcos — a H io, 

^ M H 2M 1 



(12) 
(13) 



The coherence between adjacent sites is achieved when 
w <C tN/M, which corresponds to fluctuations in the 
relative phase between neighboring sites much smaller 
than 1 (see the Appendix). If this condition is satis- 
fied, we can distinguish two limits, depending on how 
the kinetic energy per particle and the interaction energy 
per particle compare to each other: t 3> wN/M defines 
the Rabi regime, and t <C wN/M defines the Josephson 
regime. These limits were initiall y in troduced in con- 
nection with the two- well problem [l£j, where they have 
been shown to correspond to number-phase coherent and 
squeezed ground states but their generalization to lattices 
is straightforward (see 0, ^| and the Appendix). Ex- 
perimentally, achieving the Rabi regime and at the same 
time preserving the validity of the tight-binding approx- 
imation can be done by tuning the scattering length in a 
magnetic field (Feshbach resonance) rather than decreas- 
ing the depth of the potential wells. 



which is precisely a circular current state with circulation 
quantized by q. 

In obtaining the Gross-Pitaevskii equation we have as- 
sumed that there is already a phase relation established 
between the M condensates, in other words the state 
is superfluid and not fragmented (Mott insulator) 



Let us now turn to the problem of the excitation spec- 
trum. To derive the Bogoliubov - de Gennes equations 
we start with the time-dependent Gross-Pitaevskii equa- 
tion and linearize it for small fluctuations of the order 
parameter around the macroscopically occupied mode. 
We obtain 



J 



E k u k {\) = -t [u k {\ + l) + Ufc (A - 1)] + + 2Nw\ Xq {\)\ 2 ] u k (X) + Nw X 2 q (X)v k (X), 
-E k v k {\) = -t MA + 1) + v k (X - 1)] + [-ii, + 2Nw\ Xq (X)\ 2 ] v k (X) + Nw X f(X)u k (X). 

I 



(14) 
(15) 



The solution of this system is found as 
1 



Ufc(A) 
v k (X) 



M 
1 



e i(2n/M)(k+q) u ^ 



>M 



e i(2-K/M)(k-q) Vk ^ 



(16) 
(17) 



with normalization \u k \ 2 — \v k \ 2 = 1; the corresponding 
equations for u k and v k become 



2t cos — [k 




27r 

+ 2t cos —q + 

M 


N 

— w 

M 


u k 








N 

+ M WVk 


= E k u k , 


(18) 
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-q) 


— 2t cos —o — 
M 


N 

M 


v k 








N 

-M WUk 


= E k v k . 


(19) 



Solving this system gives 



(±) 2vr, . 2tt / 2N 

E k =2tsm-ksm-q±Je k U k + —w 



(20) 



where we have used the notation 

2tt / 2tt , 

e k = 2t cos — q 1 — cos — k 



(21) 



Replacing the solution for E k in Eqs. (|18I19|I we obtain 

N 



2N \ N 
£fe T \le k | £fc + -j^-w I + — w 



'I. '< \! 1 k ( '-fe + -yj- W 



N 
— i 
M 



-wu 



(±) 

k ' 

(23) 

with the restriction imposed by the normalization con- 



,(±) 



„(±) 



M 



(±) 

k ' 

(22) 



N , 



dition \u k \ 



1 which eventually will force us 



to make a choice between the two possible solutions in- 
dexed by (±). An interesting observation concerns the 
negative-energy eigenstates; assuming that we found a 

positive eigenenergy E^ with eigenvalues ( U f±) ^ ) , 

V v k ( A ) / 

the corresponding negative-energy eigenstate is — 
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with eigenvalues 



C±)* (A) 



But —E 



(±) - c;(=F) 



(±)< 
,(±)« 



(A) 
(A) 



E y \> and 



5 (A) 



where we have used 



,(±)* 



and ujj. - 1 * 



,(T). 
-fe ' 



these 



last two relations can be proved readily from Eqs. i|T8|) 
and l]19p. They show that if, say, (+) for a certain k 
is a positive-energy eigenstate, then its corresponding 
negative-energy eigenstate is a positive-energy eigenstate 
of (— ) for —k. In other words, the state with E^ rep- 
resents the "antiparticle" (in the spirit of Dirac's theory) 
of the state with E 



(+) 
-k ■ 



III. STABILITY ANALYSIS 

The stability of the states can be checked out by con- 
sidering small perturbations of the condensate state - by 
small alterations of the phase and the number of parti- 
cles on each site. It is useful to distinguish between two 
types of stabilities. 



A. Dynamical stability 

If the eigenvalues E k are complex, then the system is in 
a dynamically unstable state (or in a point of Lyapunov 
instability). The reason is that any perturbation will be 
exponentially magnified - the system tends to go as fast 
as it can as far as possible from that point. Dynamically 
stable states are those for which the condition 



2N 

e k [e k + —w 



> 



(24) 



is satisfied for any k. Let us first note that cos 2-nqjM = 
does not yield properly normalized solutions, according 
to Eqs. |22) and (J23J), so we will exclude from the begin- 
ning the states q — ±M/4. We distinguish two cases: 

1°] £fc > which is equivalent to cos (2irq/M) > 0, or 
q is in the interval (-M/4, M/4). 

2°] e^. < and — > 2wN/M. Since the minimum 
value of — £fe is reached when k = ±1, it follows that 

2tt Nw ,„ , 

cos — q< ; ; — -. (25) 

M - Mt(l -cos2tt/M) v ' 

This inequality implies also wN/M < t (1 — cos2-7r/M). 

In conclusion, the system is in a dynamically stable 
state if either condition 1°] or 2°] is satisfied. In the 
Josephson regime, only condition 1°] can be satisfied, so 
only the modes from —M/4 to M/4 are dynamically sta- 
ble. In the Rabi regime, Nw/Mt C 1 so we can distin- 
guish two cases: when M is of the order of unity, clearly 
all the modes are stable; however, when M is large, of 



the order of t/Nw, there can be dynamically unstable 
modes located at \q\ > M/4; the relative number of these 
modes, as a fraction of the total number M of mod es, is 
extremely small, since it is limited by yj Nw/Mt <C 1. 
We conclude that, in the Rabi regime, most of the states 
are dynamically stable. The dynamically unstable states 
can be treated in a linearized theory only for periods of 
time that are logarithmic with respect to the initial state 
[r?| , since the quantum fluctuations will trigger an expo- 
nentially divergent evolution away from the initial state. 



B. Thermodynamical stability 

A point of thermodynamical instability (or energetic 
instability) is a circular flow state which is still not a 
local minimum of the energy functional; however, small 
perturbations do not dynamically bring the system far 
from the initial state in the absence of dissipation. If dis- 
sipation is introduced (or, as in our case, if the system 
suffers thermalization by collisions), the system will not 
stay arbitrarily close to the initial state but instead decay 
to a thermodynamically stable state. The name "thermo- 
dynamical instability" is thus justified by the fact that 
such a state cannot be in thermodynamic equilibrium. 
The thermodynamically unstable states have a real exci- 
tation spectrum which is characterized by the existence 
of eigenenergies with negative Ek (with the normalization 
corresponding to positive eigenstates, \u k \ 2 — \v k \ 2 = 1), 
and are known to produce interesting effects in optical 
lattices |l8j . For vortices in harmonic traps they are re- 
sponsible for the precession of the vortex core around the 
condensate axis hj. This occurs because the system in 
a thermodynamically unstable vortex state can reduce 
its energy by transferring particles from the condensate 
to the negative-energy modes; for an un-pinned vortex 
for example this happens by few-particle excitations to 
the core mode [2(j. Let us now study the case of ther- 
modynamical stability in situations 1°] and 2°] described 
above. 

1°] cos2irq/M > 0. In this case only E^ is a valid 
solution. Indeed, from Eq. l(2*2l . the condition \uk\ > |«k| 
can be satisfied if 



etTl ft 4 



2N 



N 
M ? 



N 

< ~T7 W : 

M 



(26) 



and clearly the solution with the lower sign cannot satisfy 
it. 

Let us show that the upper sign solution always 
satisfies this inequality. If the expression under the 
modulus in Eq. I|2t)|l is positive, the inequality is 
fulfilled trivially. If it is negative, it becomes — 
yje k (e k + 2wN/M) > -2wN/M, or in another form 

£fe/ £fc — \Afc (ft + 2wN/M) < 1, which is satisfied for 
any values of the parameters. 

The values of the Bogoliubov amplitudes u k and v k are 
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obtained from Eq. (1221 as 



Vk 



Uh 



Ve fc + 2wN/M - ygfc 
2[e fc (e fe + 2 W 7V/M)] 1/4 
^e fc + 2wjV/Af + 



2[e fe (e fc + 2 W iV/M)] 1/4 



(27) 
(28) 



We now have to impose the condition of thermody- 
namical stability: > for any fc. This condition 
is certainly satisfied when 2nk/M is in the first and sec- 
ond quadrants, but not necessarily when sin27rfc/M < 0. 

Enforcing E { k +) > for fc < yields 



2tt , n 2 2?r ' 
1 + cos — fc — 2 cos — q 
M M 



N 2tt 

< TT WC0S TT<7- (29) 



The maximum of the left hand side (LHS) is achieved 
when fc = — 1. With some algebraic manipulations the 
constraint can be put in the form 

7T 7T N 2jT 

2tsin—(2q- 1) sin — (2q + 1) < -wcos— q. (30) 

If this condition is satisfied we have a so-called persistent 
current flowing in our geometry - this current does not 
decay in time (say by thermal fluctuations) to a lower en- 
ergy state. As it should be, the ground state (obtained for 
q = 0) trivially satisfies this inequality since in this case 
the LHS turns negative. The thermodynamically stable 
states with q are local minima of the energy func- 
tional; in their case, it costs energy to create excitations. 
We expect persistent currents to be those that will form 
via symmetry-breaking mechanisms, and survive enough 
to be finally detected. 

2°] There are no persistent currents in this case, as one 
can readily check. Indeed, this time e k < —2wN/M, and 
in order to have |ufc| 2 > |«fc| 2 , we need to impose, from 
(|22J) . the condition 



0< |e fc |± 



27V 



-w 



2N 

< n w - 



(31) 



Only the lower sign ensures the validity of these inequal- 
ities. The corresponding energy will be then E 



(-) 



E 



(-) 



2tt 



2n 



= 2t sin — -fc sin — q — 



M 



M 



,,[e k + ^w). (82) 



But now it is obvious that we cannot have Ej. ' > for 
all fc's, because if q < then the modes with fc > will 
be thermodynamically unstable, and if q > then the 
modes with fc < will not be stable. Thus, the only case 
in which we can have thermodynamical stability is case 

10 ]' 

In conclusion, the number of stable states of circular 
currents depends, in general, on the number of particles 
N, on the number of condensates M, and on the ratio t/w 



(which can be controlled experimentally by varying the 
barrier potentials between the condensates). The precise 
form of this dependence is given by the condition |q| < 
M/4, combined with Eq. I|3U|I . For example, according 
to Eq. (|30|l there cannot be stable currents with M = 3 
and M = 4, M = 5, 6, 7, 8 give at most two stable circular 
currents, M = 9, 10, 11, 12 give at most four stable ones, 
etc. 

It is instructive to see what happens in the Rabi 
regime, defined as t > wN/M, and in the Josephson 
regime, where tM/N < w < tN/M. In the Joseph- 
son regime it is clear that Eq. l|30"l) is satisfied, so we 
have both thermodynamical and dynamical stability for 
cos 27r g/A/ > 0. More interesting is what happens in the 
Rabi regime. If M is of the order of unity clearly Eq. 
(|30|l cannot be fulfilled; when M becomes of the order of 
2t/Nw, there can be stable modes with momenta given 
by 2?rg/M < y / Nw/2tM + (tt/M) 2 . They are then a 
very small fraction of the total number of modes M and 
they lead to small momenta of the circulating fluid. For 
almost all practical purposes these modes can be assim- 
ilated with q = - in a real experiment they would be 
difficult to be distinguished from the ground state. The 
fact that the rest of the modes (with larger q) in the 
Rabi regime are thermodynamically unstable can be un- 
derstood physically in a simple way. For these modes, 
we can simply put w — and search for solutions with 
iik = 1 and Vk = 0. If we look at Eq. (|22p. it is clear that 
we need to take the upper sign only into account. The 
energy in this case is 



2tt . 2vr 2tt 
Ei. = 2t sin — fc sin — o + 2t cos — q 



1 



cos fc 

M 



which can also be written in the more relevant form 

2n 2n 
E k = 2tcos — o-2icos — (k + q). (33) 

This shows that Ek is the difference between the phase- 
twisting energy per particle corresponding to a circula- 
tion q and the energy of an excitation around the macro- 
scopically occupied mode. In other words, to create an 
elementary excitation fc relative to the circulation cur- 
rent, we need to transfer an atom from the condensate 
to the single-particle state k + q. This expression is pos- 
itive for all fc's only if q = - so all the flow states are 
thermodynamically unstable. Instead, we have no dy- 
namical instability. It is interesting to note that it is 
indeed the interaction that stabilizes in the end the sys- 
tem: a large number of persistent currents cannot exist 
in its absence. Thus, phase-symmetry breaking mecha- 
nisms, for instance, can effectively work only for systems 
in which there is some extra energy to compensate for the 
loss of kinetic energy due to phase fluctuations. These 
results are summarized in Table [I] In the Fock regime 
(w 3> tN/M) we cannot write Bogoliubov equations, so 
our treatment would not be valid. 

For thermodynamically stable states (persistent cur- 
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TABLE I: Stability of circular currents: the Rabi and Joseph- 
son regimes. We find that a relatively large number (com- 
pared with M) of persistent currents can exist only above a 
certain critical value of the interaction (which turns out to 
be in between the Rabi and Josephson regimes) and only for 
certain circulations (or angular momenta). 

| REGIME | RABI | JOSEPHSON 

defined by t > wN/M tM/N < w < tN/M 

dynamically stable most \q\ < M/4 

thermodynamically stable very few \q\ < M/4 



rents) the quasiparticle spectrum obtained, 




is different from that of the ground state, which can serve 
as a method to detect a circular flow state |2l|. Also, 
similar to a phenomenon called the Sagnac effect in rel- 
ativistic physics, we note that the quasiparticles k and 
— k propagate with different speeds around the loop, be- 
cause they are excitations on top of a state with a flow 
in a certain direction. We have then a condensed-matter 
equivalent of the Sagnac effect which is related to the 
lifting of the degeneracy of the excitations k, —k due to 
the rotation of the condensate. 



IV. CONCLUSIONS 

We have obtained the excitation spectrum of a circular 
flow state in a gas of weakly interacting bosons and ana- 
lyzed its dynamical and thermodynamical stability. Test- 
ing these stability predictions, by changing the tunneling, 
the interaction, or the number of particles is within the 
reach of present-day experimental technology. 
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APPENDIX: NUMBER FLUCTUATIONS 

We give here a justification of the limits that define the 
Rabi and Josephson regimes. As in the case of the double 



well [15|, these are set respectively by the condition of 
Poissonian or sub-Poissonian fluctuations of the number 
of atoms in each well. 

Let us look now at the ground state of the system. We 
will use the same notations as throughout the main part 
of the paper with a superscript to indicate that they are 
taken at q = 0. Since all sites are equivalent, it is enough 
to look at one of them, say A = 0. If A^o is the number 
of atoms condensed on the mode q = 0, then bg ~ \/No, 
and from Eqs. (10,11) one gets 

*K0) - at(o)a(o) - ^ + ^°£(S fc + &t) 

fc^O 

+h E (A-i) 

But No = N — J2kM ^I^ fc so ignoring the smallest terms 
in the expansion the number of particles field operator 
at A = can be approximated as 

»<°>^ + f EMI). (A.2) 

fe^O 

which implies n = (n(0)) = N/M . In other words, the 
splitting of the particle number operator at any point in 
the lattice (the density operator in the limit M ^ 1) into 
a mean-field value and a fluctuating field of zero average 
results naturally from the Bogoliubov theory. Now we 
can calculate the fluctuations on the ground state 

N f° 

a* = (n 2 (0)) - (h(0)) 2 * E (A-3) 

fc^O k 

Here e° = 4t{sinnk/M) 2 and E% = y/el(e% + 2wN/M). 
If the "size" M of the system is large the sum can be 
transformed into an integral and calculated analytically. 
The result is 

2 2n nr , 

<7„ ~ — arctanW , (A. 4) 

the same as obtained in Ref. ^(| for a similar Hamil- 
tonian by using a phonon approach. The parameter 
that controls the fluctuations of the density (and phase) 
is then t/nw; the Rabi regime, with Poissonian fluc- 
tuations <7„ ~ \fn, is obtained for t 3> nw, while 
the Josephson regime, with sub-Poissonian fluctuations 
u n ~ (Stn/nw) 1 ^ 4 -C y/ri is obtained in the limit t <C nw. 
The condition a n ^> 1 (or, equivalently, phase fluctua- 
tions much smaller than 1) sets the limit of validity of the 
Bogoliubov approach: tn 3> w. In the Josephson regime, 
the excitation spectrum is ~ V 8tnw sin irk/M. In the 
Rabi regime, ~ e° for momenta 2nk/M much larger 
than a critical value J nw/t <C 1 (that is, for most of 
the modes available in the system). Below this critical 
value, if there are modes available (in other words if the 
minimum momenta allowed by quantization 2n/M is of 



7 



the order of or less than \J nw/t <C 1), the spectrum is 
phonon-likc, ~ \/8tnwTrk/M. One can regard the 
Rabi regime as being essentially a noninteracting limit 
since the fraction of excitation modes that do not have a 



spectrum as given by simply setting w = is negligible 
both in the case of M ^> 1 and M of the order of unity 
(when, since yjnwjt <C 2n/M, all the modes correspond 
in fact to a noninteracting spectrum). 
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